Spatially-dependent sensitivity of superconducting meanders 
as single-photon detectors 



G. R. Berdiyorov, M. V. Milosevic, and F. M. PeetersQ 

Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 

(Dated: March 8, 2013) 

The photo-response of a thin current-carrying superconducting stripe with a 90-degree turn is 
studied within the time-dependent Ginzburg-Landau theory. We show that the photon acting near 
the inner corner (where the current density is maximal due to the current crowding [J. R. Clem 
and K. K. Berggren, Phys. Rev. B 84, 174510 (2011)]) triggers the nucleation of superconducting 
vortices at currents much smaller than the expected critical one, but does not bring the system to 
a higher resistive state and thus remains undetected. The transition to the resistive state occurs 
only when the photon hits the stripe away from the corner due to there uniform current distribution 
across the sample, and dissipation is due to the nucleation of a kinematic vortex-antivortex pair 
near the photon incidence. We propose strategies to account for this problem in the measurements. 

PACS numbers: 74.78.Na, 73.23.-b, 74.25.Fy 



Superconducting current-carrying thin-film stripes 
have recently received a revival of interest due to their 
promising application for single-photon detection [T] with 
a high maximum count rate, broadband sensitivity, fast 
response time and low dark counts (2H1]. The single- 
photon absorption event leads to the formation of a 
nonequilibrium "hotspot" with suppressed superconduc- 
tivity, the area of which grows in time, forming a normal 
belt across the strip [5] . The latter leads to redistribution 
of the current between the now resistive superconductor 
and a parallel shunt resistor, where a voltage pulse is de- 
tected, before a hotspot region cools down (on a timescale 
of few hundred picoseconds [Hj ) and the system returns to 
its initial superconducting state. Although the hotspot 
mechanism nicely explains the photon detection in the 
visible and near UV range, the detection mechanism in 
the near-infrared range is still debated (see e.g. Ref. [7] 
and references therein). Superconducting fluctuations, 
e.g., excitation of superconducting vortices [SHU], have 
been put forward as an explanation. Dissipative cross- 
ing of such vortices, which hop over the edge barrier, or 
are created due to the unbinding of thermally activated 
vortex-antivortex pairs, provides a good description of 
the experiment [TO]. This vortex-based mechanism was 
also shown to be the dominating origin for dark counts 
[TT| \12\ fl4] . which leads to decoherence in the photon 
detection process. 

To increase the efficiency of the photon counting, de- 
tectors are usually fabricated as a meandering supercon- 
ducting wire [3]. However, these structures are vulnera- 
ble to edge imperfections [3], which significantly reduce 
the photon counting rate. Moreover, the critical current 
in these systems is mostly determined by sharp inner 
corners where the supercurrent density is maximal due 
to current crowding [T3 US]. While the appearance of 
edge imperfections can be reduced by present day tech- 
nology [3], the effect of current crowding in such me- 
andering geometries still demands further investigations. 



In this work we therefore study the effect of the turns 
of a meandering superconducting stripe on the response 
to a single-photon absorption event. Counterintuitive to 
many, we reveal that the current crowding at meander- 
ing turns does not facilitate dissipative vortex crossings 
upon the photon impact. Actually, the turning corner 
is virtually insensitive to photon absorption, which must 
be accounted for in practice. 

We consider a superconducting strip (with thickness 
d <C £, A and width w <C A = 2X 2 /d, where £, A are the 
coherence length and magnetic penetration depth) with 
a 90-degree turn in the presence of a transport current 
(see Fig. [I]). For this system we solve the following gen- 
eralized time-dependent Ginzburg-Landau equation |17j 

H ,,,,, ol + *¥> + V = V - lA ^ 
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which is coupled with the equation for the electrostatic 
potential Aip = div{lm[ip* ( V — iA)^]}. Here distance 
is scaled to £(0), the vector potential A is in units of 
$o/27r£(0)> temperature T is in units of T c , time is in 
units of to = 47rA(0) 2 /p I1 c 2 (p n is the normal state re- 
sistivity), and voltage is scaled to <pa = h/2et . Using 




FIG. 1: (Color online) The model system: a superconducting 
stripe (width w, length L and thickness gNCA,£) with a 90° 
turn. The current is applied via normal-metal contacts and 
the voltage is measured at a small distance away from these 
leads. Impact of a photon is modeled by a hotspot with radius 
R. Arrows indicate the supercurrent distribution. 
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FIG. 2: (Color online) I — V characteristics of the sample of 
size L — 400 nm and width w = 110 nm (solid black curve) 
and w = 160 nm (dashed red curve) at T = 0.9T C . Panels (1- 
4) show the snapshots of the Cooper-pair density and panels 
(1', 2') the supercurrent density for the values of the applied 
current indicated in the main panel. 



p„=18.7 fificm, £(0) = 4.2 nm and A(0) = 390 nm, which 
are typical for NbN thin films [5], one obtains io~l-25 ps 
and </Jo~0.26 mV at T = 0.9T C , which will be the working 
temperature in our simulations. The coefficients u and 7 
are chosen as u = 5.79 and 7 = 10 [17 . To model the 
thermal coupling of our sample to the substrate and the 
change of the local temperature we use the heat diffusion 
equation [T5] : 



8T 

u— = <:AT+jl-r ] (T-T ), 



(2) 



where To is the bath temperature, v is the effective heat 
capacity, £ is the effective heat conductivity coefficient, 
j n is the normal current density, and 77 is the heat trans- 
fer coefficient which governs the heat removal from the 
sample. Following the approach of Ref. [TI5] we used v — 
0.05, ( = 0.05 and 77 = 2 • 10~ 3 . This allows us to treat 
the formation and expansion of photon-induced hotspots 
in our approach. We solve the above equations self- 
consistcntly on a 2D Cartesian grid using the Euler and 
multi-grid iterative procedures. We use superconducting- 
vacuum boundary conditions (V — iA)^| n = 0, Vy>|„ = 
and VT|„ = at all sample boundaries, except at the 
current contacts where we use V = 0, T = T and 
Vy| n = —j, with j being the applied current density 
in units of j = c0 o /87r 2 £(O)A(O) 2 . 

As a representative example, we consider a supercon- 
ducting stripe as shown in Fig. [TJ of length L = 400 nm 
and width w — 110 nm, where we use the parameters 
for NbN thin films (i.e., £(0) = 4.2 nm and A(0) = 390 
nm). Time averaged voltage vs. applied current charac- 
teristics of the sample is shown in Fig. [2] for two values 
of the sample width w. With increasing the external 
current, zero resistance of the sample is maintained up 



to a threshold current density j c = 0.382jo (see solid 
black curve), above which the system goes into the re- 
sistive state. This resistive state is characterized by the 
periodic nucleation of vortices near the inner corner (see 
panels 2 and 2'), where the current density is maximal 
due to the current crowding [TS1[TS] (see panels 1 and 1'). 
However, the distribution of the supercurrent is strongly 
inhomogeneous, and decays fast away from the inner cor- 
ner towards the outer one. The Lorentz force drives the 
nucleated vortex towards the outer corner of the sample 
where it leaves the sample (panel 2), but this motion is 
slow and weakly dissipative. The nucleation rate of vor- 
tices in the inner corner increases with further increasing 
the applied current, and at sufficiently large current (la- 
beled j c d) the system transits to a higher dissipative state 
with a larger voltage jump, characterized by fast-moving 
(kinematic) vortices (see panel 3) [T21 [20] ■ We point out 
that the critical current j c decreases considerably with 
increasing width of the sample, while j c d only moder- 
ately decreases (see dashed curve in Fig. [2] for w = 160 
nm). The weakly-dissipative state is still characterized 
by vortex nucleation near the corner (see panel 4), but 
starts at low current (low j c ) and occurs in a broader 
range of the applied current (j c < j < j c d)- 

In what follows, we study the response of our system 
to a single-photon absorption event. Following the effec- 
tive temperature approach [21] . we assume that the sin- 
gle photon creates a hotspot of radius R, where the local 
temperature becomes T = 2 T c (see Ref. [1% for the ef- 
fect of such instant increase of the temperature). First, 
we consider the case when the photon acts in the middle 
of the sample away from the corner. Fig. [3|a) shows the 
V(t) characteristics of the sample with L = 400 nm and 
w = 110 nm for different values of the applied current 
density j. For each value of j, we started from the state 
obtained during the current increasing regime (i.e. states 
from Fig. [2]). The photon acts on the sample over the 
time interval of At = 25 ps (the shaded area in Fig. [3ja)), 
creating a hotspot of radius i? = 15 nm (c.f. panel 1 in 
Fig. [3]). At low bias currents (j < 0.315jo) the super- 
conducting state is stable against the photon action (not 
shown here). With increasing applied current, the sys- 
tem reacts to the photon absorption event by nucleating 
a vortex-antivortex pair (see panel 2), which is subse- 
quently unbound and split towards the edges of the sam- 
ple by the current, resulting in a voltage pulse (see solid 
curve in Fig. |3ja)). The system relaxes back into its ini- 
tial state after the vortex and the antivortex have left the 
sample and no voltage signal is observed at later times. 
The amplitude and duration of the voltage signal, as well 
as the delay in the response of the system to the photon, 
all depend on the applied current value. At larger cur- 
rent (still below j c ), several vortex-antivortex pairs can 
be formed one after another leading to several voltage 
pulses per single photon absorption (see dashed curve in 
Fig. ga)). With further increasing j, the fast moving 
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FIG. 3: (Color online) Voltage vs. time characteristics of 
the sample (of size L=400 nm and w=110 nm) at T = 0.9T C 
for different applied currents. A single-photon (with pulse 
duration At — 25 ps (shaded area in (a)) and spot radius 
R = 15 nm) acts: (a) in the middle of the stripe, c.f. panel 
1, or (b) near the inner corner, c.f. panel 4. Panels 1-7 show 
snapshots of \ip\ 2 (zoomed at the part of interest) at times 
indicated on the V(t) curves. Arrows in panel 2 show the 
direction of motion of the unbinding vortex-antivortex pair. 



vortices create a normal belt across the sample (panel 3) 
bringing the sample into the resistive state. A finite av- 
erage voltage develops across this resistive region which 
now can be detected electronically (see dotted curve in 
Fig. [3^a)). In this sense, for applied currents just be- 
low j c , we confirm the recently predicted vortex- assisted 
single-photon counting mechanism 112] . where each 
photon is detected thanks to the generated periodic mo- 
tion of kinematic vortex-antivortex pairs [13 . 

Fig. [3^b) shows the voltage response of the system 
when the photon acts near the inner corner of the sample 
(c.f. panel 4). The photoresponse of the sample becomes 
totally different from the case discussed in Fig. [3ja): (i) 
the photon generates a vortex near the inner corner of the 
sample (panel 5) at applied current values much smaller 
than the critical one (i.e., starting from j — 0.25jo). (ii) 
The amplitude/duration of the voltage signal due to the 
motion of this vortex is much smaller/longer as com- 
pared to the signal due to the vortex-antivortex unbind- 
ing (compare solid curves in Figs. |3ja) and (b)), which 
indicates the slow motion of the vortex, (iii) Instabili- 
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FIG. 4: (Color online) Phase diagram: voltage response of 
the sample of Fig. [3] as a function of the location of the 
photon absorption for two values of the applied current: (a) 
j = 0.373 jo and (b) j = 0.383j . 



ties can occur during the photon absorption (see dashed 
curve) due to the formation and dissociation of multi- 
quanta vortices (panel 6) near the inner corner, where 
supercurrent density is maximal (see panel 1' in Fig. |2j. 
(iv) No transition into the phase-slip state (similar to 
the one in inset 3 of Fig. [2]) occurs and the resistive 
state is characterized by a weakly dissipative ("cold") 
vortex crossing [TT] (see the dotted curve and inset 7 in 
Fig. [3](b)) . Therefore, although the photon acting near 
the corner does trigger the nucleation of superconduct- 
ing vortices at currents much smaller than j c , it does not 
bring the system to a higher resistive state that can be 
detected electronically. 

Our findings are summarized in Fig. [4] where we plot a 
phase diagram, showing the voltage response of the sam- 
ple as a function of the location of the photon absorption. 
The voltage is averaged over the time interval t — 0.4 — 2 
ns to exclude the oscillating features of transient voltages 
during the photon action (such as those shown in Fig. |3| . 
It is clear from this figure that the photoresponse of the 
meandering detector is spatially dependent with maximal 
sensitivity away from the turning point. Only at currents 
above j c (Fig. [4](b) ) the regions near the turning corner 
of the sample become responsive to the incident photon. 
Here we assume that the small and short-lived resistance 
due to the slow moving vortices nucleated near the inner 
corner of the sample (see solid curve in Figs. [2]and|3][b)) 
is not sufficient to induce current at the shunting load, 
so that it is beyond the sensitivity of the photodetecting 
measurement. 

In conclusion, we confirmed that vortex- assisted en- 
ergy dissipation is the dominant origin for single pho- 
ton counting in superconducting single-photon detectors. 
However, sensitivity of the commonly used meandering 
detectors is spatially dependent with maximal response 
to single-photon absorption away from the sharp corners, 
and with small response at the turns. Our suggestion is 
therefore to perform measurements at currents just under 
jc the critical current which by itself induces a moving 
vortex at the inner corners of the meander, and calcu- 
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late the photon density as the number of detected counts 
over the area of only the straight parts of the detector. 
Alternatively, one can perform measurements at currents 
just under j^, which corresponds to the jump to large 
dissipation. In this weakly dissipative regime every im- 
pact of a photon would lead to a count, but one must 
cope with enduring heating issues and non-zero dynamic 
resistance of the detector in the electronic circuit. How- 
ever, the latter strategy may be the only one suitable for 
wider meandering detectors, which provide larger surface 
to capture photons, but suffer from particularly low j c . 
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